Integrated approach to generate artificial samples with low tumor fraction for somatic variant calling benchmarking.

BACKGROUND: High-throughput sequencing (HTS) has become the gold standard approach for variant analysis in cancer research. However, somatic variants may occur at low fractions due to contamination from normal cells or tumor heterogeneity; this poses a significant challenge for standard HTS analysis pipelines. The problem is exacerbated in scenarios with minimal tumor DNA, such as circulating tumor DNA in plasma. Assessing sensitivity and detection of HTS approaches in such cases is paramount, but time-consuming and expensive: specialized experimental protocols and a sufficient quantity of samples are required for processing and analysis. To overcome these limitations, we propose a new computational approach specifically designed for the generation of artificial datasets suitable for this task, simulating ultra-deep targeted sequencing data with low-fraction variants and demonstrating their effectiveness in benchmarking low-fraction variant calling. RESULTS: Our approach enables the generation of artificial raw reads that mimic real data without relying on pre-existing data by using NEAT, a fine-grained read simulator that generates artificial datasets using models learned from multiple different datasets. Then, it incorporates low-fraction variants to simulate somatic mutations in samples with minimal tumor DNA content. To prove the suitability of the created artificial datasets for low-fraction variant calling benchmarking, we used them as ground truth to evaluate the performance of widely-used variant calling algorithms: they allowed us to define tuned parameter values of major variant callers, considerably improving their detection of very low-fraction variants. CONCLUSIONS: Our findings highlight both the pivotal role of our approach in creating adequate artificial datasets with low tumor fraction, facilitating rapid prototyping and benchmarking of algorithms for such dataset type, as well as the important need of advancing low-fraction variant calling techniques.

ARGNet: using deep neural networks for robust identification and classification of antibiotic resistance genes from sequences.

BACKGROUND: Emergence of antibiotic resistance in bacteria is an important threat to global health. Antibiotic resistance genes (ARGs) are some of the key components to define bacterial resistance and their spread in different environments. Identification of ARGs, particularly from high-throughput sequencing data of the specimens, is the state-of-the-art method for comprehensively monitoring their spread and evolution. Current computational methods to identify ARGs mainly rely on alignment-based sequence similarities with known ARGs. Such approaches are limited by choice of reference databases and may potentially miss novel ARGs. The similarity thresholds are usually simple and could not accommodate variations across different gene families and regions. It is also difficult to scale up when sequence data are increasing. RESULTS: In this study, we developed ARGNet, a deep neural network that incorporates an unsupervised learning autoencoder model to identify ARGs and a multiclass classification convolutional neural network to classify ARGs that do not depend on sequence alignment. This approach enables a more efficient discovery of both known and novel ARGs. ARGNet accepts both amino acid and nucleotide sequences of variable lengths, from partial (30-50 aa; 100-150 nt) sequences to full-length protein or genes, allowing its application in both target sequencing and metagenomic sequencing. Our performance evaluation showed that ARGNet outperformed other deep learning models including DeepARG and HMD-ARG in most of the application scenarios especially quasi-negative test and the analysis of prediction consistency with phylogenetic tree. ARGNet has a reduced inference runtime by up to 57% relative to DeepARG. CONCLUSIONS: ARGNet is flexible, efficient, and accurate at predicting a broad range of ARGs from the sequencing data. ARGNet is freely available at https://github.com/id-bioinfo/ARGNet , with an online service provided at https://ARGNet.hku.hk . Video Abstract.

Application value of metagenomic next-generation sequencing in hematological patients with high-risk febrile neutropenia.

Background: Metagenomic next-generation sequencing (mNGS) is a novel non-invasive and comprehensive technique for etiological diagnosis of infectious diseases. However, its practical significance has been seldom reported in the context of hematological patients with high-risk febrile neutropenia, a unique patient group characterized by neutropenia and compromised immune responses. Methods: This retrospective study evaluated the results of plasma cfDNA sequencing in 164 hematological patients with high-risk febrile neutropenia. We assessed the diagnostic efficacy and clinical impact of mNGS, comparing it with conventional microbiological tests. Results: mNGS identified 68 different pathogens in 111 patients, whereas conventional methods detected only 17 pathogen types in 36 patients. mNGS exhibited a significantly higher positive detection rate than conventional methods (67.7% vs. 22.0%, P < 0.001). This improvement was consistent across bacterial (30.5% vs. 9.1%), fungal (19.5% vs. 4.3%), and viral (37.2% vs. 9.1%) infections (P < 0.001 for all comparisons). The anti-infective treatment strategies were adjusted for 51.2% (84/164) of the patients based on the mNGS results. Conclusions: mNGS of plasma cfDNA offers substantial promise for the early detection of pathogens and the timely optimization of anti-infective therapies in hematological patients with high-risk febrile neutropenia.

Pragmatic nationwide master observational trial based on genomic alterations in advanced solid tumors: KOrean Precision Medicine Networking Group Study of MOlecular profiling guided therapy based on genomic alterations in advanced Solid tumors (KOSMOS)-II study protocol KCSG AL-22-09.

BACKGROUND: Next-generation sequencing (NGS) has been introduced to many Korean institutions to support molecular diagnostics in cancer since 2017, when it became eligible for reimbursement by the National Health Insurance Service. However, the uptake of molecularly guided treatment (MGT) based on NGS results has been limited because of stringent regulations regarding prescriptions outside of approved indications, a lack of clinical trial opportunities, and limited access to molecular tumor boards (MTB) at most institutions. The KOSMOS-II study was designed to demonstrate the feasibility and effectiveness of MGT, informed by MTBs, using a nationwide precision medicine platform. METHODS: The KOSMOS-II trial is a large-scale nationwide master observational study. It involves a framework for screening patients with metastatic solid tumors for actionable genetic alterations based on local NGS testing. It recommends MGT through a remote and centralized MTB meeting held biweekly. MGT can include one of the following options: Tier 1, the therapeutic use of investigational drugs targeting genetic alterations such as ALK, EGFR, ERBB2, BRAF, FH, ROS1, and RET, or those with high tumor mutational burden; Tier 2, comprising drugs with approved indications or those permitted for treatment outside of the indications approved by the Health Insurance Review and Assessment Service of Korea; Tier 3, involving clinical trials matching the genetic alterations recommended by the MTB. Given the anticipated proportion of patients receiving MGT in the range of 50% ± 3.25%, this study aims to enroll 1,000 patients. Patients must have progressed to one or more lines of therapy and undergone NGS before enrollment. DISCUSSION: This pragmatic master protocol provides a mass-screening platform for rare genetic alterations and high-quality real-world data. Collateral clinical trials, translational studies, and clinico-genomic databases will contribute to generating evidence for drug repositioning and the development of new biomarkers. TRIAL REGISTRATION: NCT05525858.

Identification of the novel HLA-DQA1*02:32 allele by next-generation sequencing.

HLA-DQA1*02:32 differs from DQA1*02:01:01 by one nucleotide substitution in codon 210 in exon 4.

Characterisation of the novel HLA-C*15:274 allele using short and long-read sequencing technologies.

The novel HLA-C*15:274 allele, first described in a potential bone marrow donor from Brazil.

Multi-centre analytical performance verification of an IVD assay to quantify donor-derived cell-free DNA in solid organ transplant recipients.

Donor-derived cell-free DNA (dd-cfDNA) has been widely studied as biomarker for non-invasive allograft rejection monitoring. Earlier rejection detection enables more prompt diagnosis and intervention, ultimately improving patient treatment and outcomes. This multi-centre study aims to verify analytical performance of a next-generation sequencing-based dd-cfDNA assay at end-user environments. Three independent laboratories received the same experimental design and 16 blinded samples to perform cfDNA extraction and the dd-cfDNA assay workflow. dd-cfDNA results were compared between sites and against manufacturer validation to evaluate concordance, reproducibility, repeatability and verify analytical performance. A total of 247 sample libraries were generated across 18 runs, with completion time of <24 h. A 96.0% first pass rate highlighted minimal failures. Overall observed versus expected dd-cfDNA results demonstrated good concordance and a strong positive correlation with linear least squares regression r2 = 0.9989, and high repeatability and reproducibility within and between sites, respectively (p > 0.05). Manufacturer validation established limit of blank 0.18%, limit of detection 0.23% and limit of quantification 0.23%, and results from independent sites verified those limits. Parallel analyses illustrated no significant difference (p = 0.951) between dd-cfDNA results with or without recipient genotype. The dd-cfDNA assay evaluated here has been verified as a reliable method for efficient, reproducible dd-cfDNA quantification in plasma from solid organ transplant recipients without requiring genotyping. Implementation of onsite dd-cfDNA testing at clinical laboratories could facilitate earlier detection of allograft injury, bearing great potential for patient care.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data.

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

CAREx: context-aware read extension of paired-end sequencing data.

BACKGROUND: Commonly used next generation sequencing machines typically produce large amounts of short reads of a few hundred base-pairs in length. However, many downstream applications would generally benefit from longer reads. RESULTS: We present CAREx-an algorithm for the generation of pseudo-long reads from paired-end short-read Illumina data based on the concept of repeatedly computing multiple-sequence-alignments to extend a read until its partner is found. Our performance evaluation on both simulated data and real data shows that CAREx is able to connect significantly more read pairs (up to 99 % for simulated data) and to produce more error-free pseudo-long reads than previous approaches. When used prior to assembly it can achieve superior de novo assembly results. Furthermore, the GPU-accelerated version of CAREx exhibits the fastest execution times among all tested tools. CONCLUSION: CAREx is a new MSA-based algorithm and software for producing pseudo-long reads from paired-end short read data. It outperforms other state-of-the-art programs in terms of (i) percentage of connected read pairs, (ii) reduction of error rates of filled gaps, (iii) runtime, and (iv) downstream analysis using de novo assembly. CAREx is open-source software written in C++ (CPU version) and in CUDA/C++ (GPU version). It is licensed under GPLv3 and can be downloaded at ( https://github.com/fkallen/CAREx ).

Development, validation and application of single molecule molecular inversion probe based novel integrated genetic screening method for 29 common lysosomal storage disorders in India.

BACKGROUND: Current clinical diagnosis pathway for lysosomal storage disorders (LSDs) involves sequential biochemical enzymatic tests followed by DNA sequencing, which is iterative, has low diagnostic yield and is costly due to overlapping clinical presentations. Here, we describe a novel low-cost and high-throughput sequencing assay using single-molecule molecular inversion probes (smMIPs) to screen for causative single nucleotide variants (SNVs) and copy number variants (CNVs) in genes associated with 29 common LSDs in India. RESULTS: 903 smMIPs were designed to target exon and exon-intron boundaries of targeted genes (n = 23; 53.7 kb of the human genome) and were equimolarly pooled to create a sequencing library. After extensive validation in a cohort of 50 patients, we screened 300 patients with either biochemical diagnosis (n = 187) or clinical suspicion (n = 113) of LSDs. A diagnostic yield of 83.4% was observed in patients with prior biochemical diagnosis of LSD. Furthermore, diagnostic yield of 73.9% (n = 54/73) was observed in patients with high clinical suspicion of LSD in contrast with 2.4% (n = 1/40) in patients with low clinical suspicion of LSD. In addition to detecting SNVs, the assay could detect single and multi-exon copy number variants with high confidence. Critically, Niemann-Pick disease type C and neuronal ceroid lipofuscinosis-6 diseases for which biochemical testing is unavailable, could be diagnosed using our assay. Lastly, we observed a non-inferior performance of the assay in DNA extracted from dried blood spots in comparison with whole blood. CONCLUSION: We developed a flexible and scalable assay to reliably detect genetic causes of 29 common LSDs in India. The assay consolidates the detection of multiple variant types in multiple sample types while having improved diagnostic yield at same or lower cost compared to current clinical paradigm.

Direct Analysis of HIV mRNA m6A Methylation by Nanopore Sequencing.

The post-transcriptional processing and chemical modification of HIV RNA are understudied aspects of HIV virology, primarily due to the limited ability to accurately map and quantify RNA modifications. Modification-specific antibodies or modification-sensitive endonucleases coupled with short-read RNA sequencing technologies have allowed for low-resolution or limited mapping of important regulatory modifications of HIV RNA such as N6-methyladenosine (m6A). However, a high-resolution map of where these sites occur on HIV transcripts is needed for detailed mechanistic understanding. This has recently become possible with new sequencing technologies. Here, we describe the direct RNA sequencing of HIV transcripts using an Oxford Nanopore Technologies sequencer and the use of this technique to map m6A at near single nucleotide resolution. This technology also provides the ability to identify splice variants with long RNA reads and thus, can provide high-resolution RNA modification maps that distinguish between overlapping splice variants. The protocols outlined here for m6A also provide a powerful paradigm for studying any other RNA modifications that can be detected on the nanopore platform.

Dual RNA-seq Analysis of Patients' Cells and Viral Genome After Measles Infection.

During the infection of a host cell by an infectious agent, a series of gene expression changes occurs as a consequence of host-pathogen interactions. Unraveling this complex interplay is the key for understanding of microbial virulence and host response pathways, thus providing the basis for new molecular insights into the mechanisms of pathogenesis and the corresponding immune response. Dual RNA sequencing (dual RNA-seq) has been developed to simultaneously determine pathogen and host transcriptomes enabling both differential and coexpression analyses between the two partners as well as genome characterization in the case of RNA viruses. Here, we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on - but not restricted to - measles virus (MeV) as a pathogen of interest. The application of dual RNA-seq technologies in MeV-infected patients can potentially provide valuable information on the structure of the viral RNA genome and on cellular innate immune responses and drive the discovery of new targets for antiviral therapy.

Enzyme-assisted high throughput sequencing of an expanded genetic alphabet at single base resolution.

With just four building blocks, low sequence information density, few functional groups, poor control over folding, and difficulties in forming compact folds, natural DNA and RNA have been disappointing platforms from which to evolve receptors, ligands, and catalysts. Accordingly, synthetic biology has created "artificially expanded genetic information systems" (AEGIS) to add nucleotides, functionality, and information density. With the expected improvements seen in AegisBodies and AegisZymes, the task for synthetic biologists shifts to developing for expanded DNA the same analytical tools available to natural DNA. Here we report one of these, an enzyme-assisted sequencing of expanded genetic alphabet (ESEGA) method to sequence six-letter AEGIS DNA. We show how ESEGA analyses this DNA at single base resolution, and applies it to optimized conditions for six-nucleotide PCR, assessing the fidelity of various DNA polymerases, and extending this to AEGIS components with functional groups. This supports the renewed exploitation of expanded DNA alphabets in biotechnology.

Hybracter: enabling scalable, automated, complete and accurate bacterial genome assemblies.

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants. They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.

Identification of five new HLA-C alleles by next-generation sequencing.

Five novel HLA-C alleles detected by next-generation sequencing: HLA-C*02:02:73, -C*03:04:106, -C*06:382, -C*07:1114Q and -C*12:408.

Characterisation of the novel HLA-B*35:593 allele using short and long-read sequencing technologies.

The novel HLA-B*35:593 allele, first described in a potential bone marrow donor from Brazil.

Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing.

Whole-genome reconstruction of bacterial pathogens has become an important tool for tracking transmission and antimicrobial resistance gene spread, but highly accurate and complete assemblies have largely only historically been achievable using hybrid long- and short-read sequencing. We previously found the Oxford Nanopore Technologies (ONT) R10.4/kit12 flowcell/chemistry produced improved assemblies over the R9.4.1/kit10 combination, however long-read only assemblies contained more errors compared to Illumina-ONT hybrid assemblies. ONT have since released an R10.4.1/kit14 flowcell/chemistry upgrade and recommended the use of Bovine Serum Albumin (BSA) during library preparation, both of which reportedly increase accuracy and yield. They have also released updated basecallers trained using native bacterial DNA containing methylation sites intended to fix systematic basecalling errors, including common adenosine (A) to guanine (G) and cytosine (C) to thymine (T) substitutions. To evaluate these improvements, we successfully sequenced four bacterial reference strains, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, and nine genetically diverse E. coli bloodstream infection-associated isolates from different phylogroups and sequence types, both with and without BSA. These sequences were de novo assembled and compared against Illumina-corrected reference genomes. In this small evaluation of 13 isolates we found that nanopore long-read-only R10.4.1/kit 14 assemblies with updated basecallers trained using bacterial methylated DNA produce accurate assemblies with ≥40×depth, sufficient to be cost-effective compared with hybrid ONT/Illumina sequencing in our setting.

Genomic profiling of pleomorphic rhabdomyosarcoma reveals a genomic signature distinct from that of embryonal rhabdomyosarcoma.

Pleomorphic rhabdomyosarcoma (PRMS) is a rare and highly aggressive sarcoma, occurring mostly in the deep soft tissues of middle-aged adults and showing a variable degree of skeletal muscle differentiation. The diagnosis is challenging as pathologic features overlap with embryonal rhabdomyosarcoma (ERMS), malignant Triton tumor, and other pleomorphic sarcomas. As recurrent genetic alterations underlying PRMS have not been described to date, ancillary molecular diagnostic testing is not useful in subclassification. Herein, we perform genomic profiling of a well-characterized cohort of 14 PRMS, compared to a control group of 23 ERMS and other pleomorphic sarcomas (undifferentiated pleomorphic sarcoma and pleomorphic liposarcoma) using clinically validated DNA-targeted Next generation sequencing (NGS) panels (MSK-IMPACT). The PRMS cohort included eight males and six females, with a median age of 53 years (range 31-76 years). Despite similar tumor mutation burdens, the genomic landscape of PRMS, with a high frequency of TP53 (79%) and RB1 (43%) alterations, stood in stark contrast to ERMS, with 4% and 0%, respectively. CDKN2A deletions were more common in PRMS (43%), compared to ERMS (13%). In contrast, ERMS harbored somatic driver mutations in the RAS pathway and loss of function mutations in BCOR, which were absent in PRMS. Copy number variations in PRMS showed multiple chromosomal arm-level changes, most commonly gains of chr17p and chr22q and loss of chr6q. Notably, gain of chr8, commonly seen in ERMS (61%) was conspicuously absent in PRMS. The genomic profiles of other pleomorphic sarcomas were overall analogous to PRMS, showing shared alterations in TP53, RB1, and CDKN2A. Overall survival and progression-free survival of PRMS were significantly worse (p < 0.0005) than that of ERMS. Our findings revealed that the molecular landscape of PRMS aligns with other adult pleomorphic sarcomas and is distinct from that of ERMS. Thus, NGS assays may be applied in select challenging cases toward a refined classification. Finally, our data corroborate the inclusion of PRMS in the therapeutic bracket of pleomorphic sarcomas, given that their clinical outcomes are comparable.

HPV genotyping in clinical samples using long-read single-molecule real-time sequencing.

Human papillomavirus (HPV) genotyping is widely used, particularly in combination with high-risk (HR) HPV tests for cervical cancer screening. We developed a genotyping method using sequences of approximately 800 bp in the E6/E7 region obtained by PacBio single molecule real-time sequencing (SMRT) and evaluated its performance against MY09-11 L1 sequencing and after the APTIMA HPV genotyping assay. The levels of concordance of PacBio E6/E7 SMRT sequencing with MY09-11 L1 sequencing and APTIMA HPV genotyping were 100% and 90.8%, respectively. The sensitivity of PacBio E6/EA7 SMRT was slightly greater than that of L1 sequencing and, as expected, lower than that of HR-HPV tests. In the context of cervical cancer screening, PacBio E6/E7 SMRT is then best used after a positive HPV test. PacBio E6/E7 SMRT genotyping is an attractive alternative for HR and LR-HPV genotyping of clinical samples. PacBio SMRT sequencing provides unbiased genotyping and can detect multiple HPV infections and haplotypes within a genotype.